Development of a simple polymer-based sensor for detection of the Pirimicarb pesticide

In this study, a sensitive and selective fluorescent chemosensor was developed for the determination of pirimicarb pesticide by adopting the surface molecular imprinting approach. The magnetic molecularly imprinted polymer (MIP) nanocomposite was prepared using pirimicarb as the template molecule, CuFe2O4 nanoparticles, and graphene quantum dots as a fluorophore (MIP-CuFe2O4/GQDs). It was then characterized using X-ray diffraction (XRD) technique, Fourier transforms infrared (FT-IR) spectroscopy, scanning electron microscope (SEM), and transmission electron microscopy (TEM). The response surface methodology (RSM) was also employed to optimize and estimate the effective parameters of pirimicarb adsorption by this polymer. According to the experimental results, the average particle size and imprinting factor (IF) of this polymer are 53.61 nm and 2.48, respectively. Moreover, this polymer has an excellent ability to adsorb pirimicarb with a removal percentage of 99.92 at pH = 7.54, initial pirimicarb concentration = 10.17 mg/L, polymer dosage = 840 mg/L, and contact time = 6.15 min. The detection of pirimicarb was performed by fluorescence spectroscopy at a concentration range of 0–50 mg/L, and a sensitivity of 15.808 a.u/mg and a limit of detection of 1.79 mg/L were obtained. Real samples with RSD less than 2 were measured using this chemosensor. Besides, the proposed chemosensor demonstrated remarkable selectivity by checking some other insecticides with similar and different molecular structures to pirimicarb, such as diazinon, deltamethrin, and chlorpyrifos.


Apparatus
The characteristics and surface properties of the pirimicarb molecularly imprinted polymer (MIP-CuFe 2 O 4 / GQDs) have been evaluated using various analyses.In this study, Fourier transform infrared (FT-IR) spectroscopy (Tensor 27, Bruker) was performed to identify the functional groups and their interactions with pirimicarb pesticide.The crystalline structures of prepared materials were analyzed by X-ray diffraction (XRD) technique (Tongda, model TD-3700, china).Moreover, the scanning electron microscope (SEM) was utilized to investigate the surface topography and composition of the aforementioned nanocomposite (TESCAN MIRA3 FEG-SEM).Transmission Electron Microscopy (TEM) images of CuFe 2 O 4 and MIP-CuFe 2 O 4 /GQDs was recorded using Philips EM 208S.The residual pirimicarb concentrations at any step were determined with a UV-VIS spectrophotometer (Specord 250, Analytik Jena), and a fluorescence spectrophotometer (LS45, PerkinElmer) was applied for the sensing process.

Preparation of CuFe 2 O 4
For a typical synthesis of CuFe 2 O 4 nanoparticles, initially Cu(NO 3 ) 2 (1870 mg) and Fe(NO 3 ) 3 (4820 mg) were dissolved in deionized water (100 mL) in 1:2 ratios, followed by 1 h of stirring at 90 °C temperature under the inert gas (N 2 ).Then, NaOH solution (30 ml, 5 M) was added dropwise to the resulting solution under stirring.After 1 h, the obtained precipitate was separated using an external magnet (3T) and washed several times with deionized water, and dried in the oven at 90 °C for 24 h.

Preparation of GQDs
The graphene quantum dots were synthesized through the bottom-up technique by pyrolysis of citric acid following the protocols of the previous report 24 .Briefly, 2000 mg of citric acid monohydrate was melted at 200 °C for 45 min in the oven and then added dropwise to 100 ml of NaOH solution with a 10 mg/ml concentration.Then, the resultant mixture was neutralized with 0.1 M NaOH after 10 min.

Synthesis procedure of polymers
In order to prepare the pirimicarb imprinted polymer, first 50 mg of pirimicarb was dissolved in chloroform (30 ml).Afterward, 100 mg of CuFe 2 O 4 nanoparticles and 100 mg of GQDs were dispersed in it, respectively.(10 min, by ultrasonic dispersion) Then 135 µL of MAA as a functional monomer, EGDMA as a cross-linking monomer (755 µL), (the ratio of template: functional monomer: cross linker = 1:8:20) and AIBN as an initiator (16 mg) were added to the obtained solution, respectively.Meanwhile, degassing was carried out before the addition of AIBN for 3 min and about 5 min after its addition.The final product (915 mg) was obtained after 19 h and at a temperature of 50 °C.Subsequently, the resultant magnetic polymer was extracted from the mixture by an external magnet.Then, it was washed with a solution including 90% acetone and 10% acetic acid by the Soxhlet method to remove the template molecule (pirimicarb) from the produced imprinted polymer construction.Nonimprinted polymer (NIP) was prepared (668 mg) in the same method as the MIP nanocomposite but without the addition of the pirimicarb.In this study, the polymer was prepared using MAA as the functional monomer and EDMA as the cross linking monomer.Hydrogen bonds are expected to be formed between the nitrogen atoms of pirimicarb and MAA as a key interaction for binding-site construction 3,[25][26][27] .Also, in this study, we have synthesized several polymers based on different ratio of template, functional monomer, and cross linker such as 1:4:20, 1:8:20, and 1:16:20.Finally, we found mentioned ratio (1:8:20) as the best value.

Assessment of effective adsorption parameters
To study the adsorption process of pirimicarb by MIP-CuFe 2 O 4 /GQDs nanocomposite, pirimicarb solutions of various values, including 10, 20, 30, 40, and 50 mg/L, were provided.pH was adjusted using 0.1 M NaOH or HCl (3 to 11).Then, specific dosages of MIP-CuFe 2 O 4 /GQDs (250 to 1250 mg/L) were added to the solutions, and the resulting compounds were stirred for 5 to 25 min.In all samples, after the pirimicarb adsorption mechanism, the MIP-CuFe 2 O 4 /GQDs were removed from the solution with an external magnet, and the concentration of the remaining pesticide in the solution was determined by a UV-1700 spectrophotometer at 313 nm.According to the results, the adsorption yield of pirimicarb by the nanocomposite and the adsorption capacity were calculated by applying Eqs.(1) and ( 2), respectively 28 .www.nature.com/scientificreports/In the equations above, R is the adsorption percentage, q e is the amount of adsorption at equilibrium (g/kg), and C i and C e are the initial and equilibrium concentrations of pirimicarb (mg/L), respectively.V is the liquid volume (L), and m is the amount of polymer (g).

Response surface methodology (RSM) experiments
To optimize pirimicarb adsorption by MIP-CuFe 2 O 4 /GQDs nanocomposite, a four-factor based model was designed using RSM (CCD) with Design Expert 11 software.For four independent variables (k), namely pH, time, initial pirimicarb concentration, and polymer quantity, the 28 (2 k + 2k + 4 = 28) experiments suggested by this model were carried out in this study.
Table 1 shows the ranges, levels, and design matrix of chosen parameters.In addition, the model suggested by response surface methodology utilizing analysis of variance was evaluated.This model contained predictor variables and demonstrated the correlation between independent and dependent variables.Significant and nonsignificant impacts of independent variables were investigated using an analysis of variance and a P-value.The model predictions were then compared to the experimental results and the correlation coefficient (R 2 ) was calculated.Further relevant statistics for analyzing model performance and efficiency, such as standard deviation (SD), coefficient of variation (CV %), desirability function (DF), and Adequate precision (AP), were also investigated.

Fluorescence measurements and sensitivity determination
To evaluate the fluorescence performance, several solutions were prepared at a pirimicarb concentration range of 0 to 50 mg/L (10 ml) and optimal pH.Then 1 mg of MIP-CuFe 2 O 4 /GQDs was dispersed (5 min) in all of them.The fluorescence intensity of each sample with an excitation wavelength of 370 nm was obtained in 10 s with a fluorescence spectrophotometer at 375-450 nm wavelengths.Then, the limit of detection (LOD) and limit of quantitation (LOQ) were calculated from the calibration line at low concentrations by the following equations (Eqs.( 3) and ( 4)) 29 .In order to check the ability of the sensor, real river water samples (Table S3) with the addition several concentration of pirimicarb (0 to 50 mg/L) were investigated.Then, based on the measured fluorescence intensity and the sensor response equation, the concentration of pirimicarb was measured.Relative Standard Deviation (RSD) for this measurement was calculated and then the accuracy of this method based on RSD was investigated.
In these equations, S a and b are the standard deviation of the intercept and the slope of the regression line, respectively.The response equation and linear model (approximately) of the sensor were also obtained based on the results of this section.Then, the sensitivity of the sensor was determined from the slope of the approximately linear model.

Selectivity and binding study
To investigate the specificity of the NIP-CuFe 2 O 4 /GQDs and MIP-CuFe 2 O 4 /GQDs nanocomposites, the binary solutions of three other different pesticides (namely Diazinon, Deltamethrin, and Chlorpyrifos) with pirimicarb (10 mg/L) were prepared.Then, 20 mg dosages of nanocomposites were added to the solutions separately, and after 1 h of stirring, they were removed from the solutions with an external magnet.The adsorption experiments were performed under identical experimental conditions for all six solutions, and the concentration of residual pesticides in solutions was measured by spectrophotometry.According to the results, the distribution coefficients of each pesticide were obtained using Eq. ( 5). (1) (2) www.nature.com/scientificreports/K d is the pesticide distribution proportion (L/g) between the adsorbent and the aqueous phase, q e is the quantity of adsorption at equilibrium (g/kg), and C e denotes the pesticide's equilibrium concentration (mg/L).
The selectivity coefficient of each adsorbent (K) for pirimicarb with regard to the competitive pesticides were calculated by Eq. ( 7).The proportion of K MIP to K NIP (Eq.( 7)) was then used to calculate the relative selectivity coefficient (K′) 30 .
Moreover, binding tests and Scatchard analysis were applied to assess the binding site and calculate the imprinting factor of the MIP.Pirimicarb binding experiments were carried out with a 10 mg dosage of each nanocomposite at a pirimicarb concentration range of 10 mg/L to 50 mg/L under optimal conditions.The amount of pirimicarb adsorbed on the polymers was calculated to be applied in the Scatchard equation (Eq.( 8)) through measuring the ultimate and initial concentration of pirimicarb with UV spectrophotometry.
In this equation, q is the binding/adsorption capacity of nanocomposites, and C e is the equilibrium concentration of pirimicarb.The K D and q max are defined as the equilibrium association constant and apparent capacity, which were obtained from the slope and intercept of the line graph wherein q/C e is plotted versus q, respectively 31 .
Following that, by computing the ratio of the apparent capacity of MIP to NIP, the imprinting factor (IF) was obtained 32 .

Reusability
The reusability of adsorbents is an essential and effective parameter for sensor material.The adsorption process was performed in optimal conditions, and desorption was done at pH 1 to 13 of Britton-Robinson buffer.The best solution was selected based on the desorption percentage (Calculated from Eq. ( 9)) 33 .Afterward, the adsorption test was performed.Then the adsorbent was separated from the solution and dried, and a desorption test was performed.This process was repeated for up to 15 cycles.

Characterization of nanocomposite and polymers
The FT-IR spectra of CuFe 2 O 4 nanoparticles, GQD, NIP-CuFe 2 O 4 /GQDs, and MIP-CuFe 2 O 4 /GQDs nanocomposites are shown in Fig. 1a.According to the spectrum of CuFe 2 O 4 , there are two main absorption bands at about 493 and 703 cm −1 , which are attributed to stretching vibrations of the octahedral and tetrahedral sites in spinels (Cu-O and Fe-O bonds), respectively 34 .Moreover, the two peaks at 1633 and 3416 cm −1 represent the adsorption peak of C=O and stretching and bending vibrations of O-H groups.The FT-IR spectrum of GQDs was shown in Fig. 1a, there is a broad peak at 3453 cm −1 is related to -OH and -NH groups stretching vibration.Also, the well-defined peak at 1585 cm −1 is assigned to bending vibrations of C=C group.The peaks at 1389 cm −1 , and 1077 cm −1 were related to the C−O (carboxy), and C−O (alkoxy) functional groups, respectively 35 .Besides, the spectra of the MIP-CuFe 2 O 4 /GQDs indicate three additional significant absorption peaks at 2989, 1731, and 1156 cm −1 , which respectively represent the presence of the methylene group in both MAA and EGDMA, stretching vibration of C=O and C-O bonds [36][37][38][39] .Also, According to the FT-IR spectrum of the pirimicarb, the peaks of this pesticide overlap with NIP and MIP nanocomposites.There is only one peak in the region of 450 cm −1 , which can also be seen in the polymer before removing the template from it.The peak in the region of 1600 cm −1 is also present in the spectrum related to polymer, but its intensity has been changed.Basically, in the synthesis of molecular imprinted polymer, the amount of template molecule is less and it is quite probable that the intensity of its peaks is low.In addition, the peak observed in the region of 2900 cm −1 in the MIP-CuFe 2 O 4 / GQDs-with pirimicarb is high intensive than the NIP-CuFe 2 O 4 /GQDs, which can be caused by the presence of pirimicarb in the structure of the MIP with pirimicarb.Also it has been decreased in intensity after washing in MIP-CuFe 2 O 4 /GQDs 27 .
Figure 1b presents the XRD patterns of the CuFe 2 O 4 , GQDs, and MIP-CuFe 2 O 4 /GQDs nanocomposite.As can be observed, the XRD pattern of the CuFe 2 O 4 indicates eight definite peaks in the range of 2θ, which are indexed to (111), ( 220), (311), ( 222), (400), (511), (440), and (533) crystallographic planes of CuFe 2 O 4 , respectively 40,41 .These peaks can also be seen in the pattern of the MIP-CuFe 2 O 4 /GQDs with a slight shift and discrepancy in the intensity, which confirms that the crystalline phase of CuFe 2 O 4 nanoparticles is not altered by MIP synthesis 42 .Besides, in the XRD pattern of GQDs, a broad peak appeared at around 2θ = 33.75concerning the (002) plane 43 , and its slight effect on the XRD pattern of the nanocomposite can be observed.For both CuFe 2 O 4 nanoparticles ( 5)  www.nature.com/scientificreports/

Results of RSM (CCD)
Experiments were prepared and conducted following the surface response technique and central composite design, and the results were evaluated using design expert software.Initially, the assumption of data normality was investigated using a normal distribution graph and a Box-Cox graph.According to the normal distribution diagram (Fig. S1-a)(see Supplementary Information (SI)), the data is distributed around the line and does not exhibit the "S" pattern, which confirming their assumption of normality 46 .The Box-Cox diagram (Fig. S1-b) (see SI) indicates that the best lambda value in the 95% confidence range of 0.12 to 1.49 is 0.77, and its current  value is 1.Hence, it is within the range of confidence, and the data does not require a transfer function to be normalized.The modeling results reveal that the prediction model is a quadratic equation depending on the factors investigated.The developed model was statistically evaluated to identify the optimal quadratic model parameters for experimental data as a significant model using analysis of variance (ANOVA) 47 .As a consequence, based on the data acquired from the response prediction model and its variables, the final equation in terms of coded factors (a quadratic polynomial function) is as follows (Eq.( 10)): The coded factors X 1 , X 2 , X 3 , and X 4 in the above equation refer to pH, time, initial concentration of pirimicarb, and nanocomposite dosage, respectively.The ANOVA analysis results of the obtained quadratic model (Eq.( 10)) are presented in Table 2.A significant model should have a non-significant lack-of-fit (p-value > 0.05) and a p-value less than 0.05 48,49 .A large F value and a low P value suggest the great significance of the related parameters 50 .As can be seen in Table 2, the F-value of 63.40 and corresponding p-value less than 0.0001 for the quadratic model confirm its significance.The non-significant lack-of-fit term (F-value = 2.63 and p-value > 0.05), implying that the residual error is not significant.Furthermore, summary statistics show R 2 , average value, SD, and CV% are equal to 0.96, 71.68, 3.87, and 5.40%, respectively.The obtained model is reproducible with regard to the coefficient of variation, which is less than 10%.Also, the predicted R 2 of 0.90 and adjusted R 2 of 0.95 are in good agreement, which demonstrate the model's incredible ability to predict experimental results.The AP statistic value was estimated to be 30.77,which represents the noise signal rate.About this value, a number higher than 4 indicates the model's appropriateness 51,52 .
Equation (10) displays the coefficients of all effective parameters.Since this equation does not provide enough information on the effect of each variable and its operation, the perturbation graph of independent variables was investigated (Fig. S2-a) (see SI).This graph indicates the impact of significant variables by demonstrating variations in the response of each factor as it moves away from the reference point (zero coded level of each factor), while all other factors remain constant at the reference point 53 .In the curve, the effect of the squared terms for all variables, except adsorbent dosage, can be observed.In addition, the slope of the graph related to time is lower in comparison to other factors, indicating that the response is less sensitive to this variable.In contrast, other variables have relatively similar slopes.Besides, increasing the initial concentration of pirimicarb decreases the adsorption process's effectiveness, whereas raising the other three influential factors leads to an increase in adsorption efficiency.Because, raising the initial concentration of pirimicarb causes its molecules to compete for placement on the polymer surface, which reduces adsorption efficiency 30 .On the Pareto chart (Fig. S2-b) (see SI), in addition to the importance of the effects of significant variables, their interactions can also be seen.As regards this chart, the polymer dosage (X 4 ) has the most significant percentage of effectiveness, and time squared has the least percentage (X 2 2 ) in adsorption efficiency.The three-dimensional graphs of the response surface of the adsorbed quantity of pirimicarb relative to the four factors were employed to get knowledge of variable interactions and evaluate the optimum amount of each variable to achieve maximum adsorption.Figure (S3-a) (see SI) depicts surface plots of pirimicarb adsorption percentage versus the pH and the polymer dosage by keeping contact time constant at 15 min and the initial concentration at 30 mg/L.In this graph, the maximum efficiency associated with the influence of two variables (pH and the polymer dosage) hit a peak at a pH of approximately 7 and the highest level of nanocomposite dosage (1250 mg/L).In Fig. (S3-b) (see SI), which demonstrates interactions between pH and the initial pirimicarb ( www.nature.com/scientificreports/concentration in a constant time (15 min) and polymer dosage (750 mg/L), the maximum efficiency associated with the influence of these two variables reached its peak at a concentration of around 10 mg/L and the highest level of pH 7 to 11.The parameters were optimized to obtain the conditions wherein all variables are ideal, and the adsorption percentage is highest.Figure S4 (see SI) depicts graphs of the actual quantity of each independent variable in which the optimal points are indicated, as well as graphs of the maximum and minimum adsorption percentage.The best parameters for removing 100 percent of pirimicarb are an initial pirimicarb concentration of 10.17 mg/L, a pH of 7.54, a contact time of 6.15 min, and a polymer dosage of 840 mg/L.This optimum condition was investigated experimentally, where 99.92 adsorption percent was attained, which has a slight discrepancy in comparison with the model estimated quantity.It can be inferred that the central composite design can efficiently characterize and estimate the optimal adsorption process circumstances for the pirimicarb by MIP-CuFe 2 O 4 / GQDs nanocomposite.

Fluorescence performance and detection sensitivity and linear range of sensor
To evaluate the detection efficiency of the proposed sensor, specific dosages of the MIP-CuFe 2 O 4 /GQDs were separately dispersed in different concentrations of pirimicarb.Then the fluorescence emission spectra of these samples were obtained under the optimal experimental condition.Also, the river water samples, whose specifications are given in Table S3, was analyzed using this method, and the results are presented in Fig. 4b.For the evaluation of the sensor's detection sensitivity, the correlation of fluorescence intensity with pirimicarb concentration was conducted, and the results are provided in Fig. 4. As can be observed, the fluorescence signal gradually increase as the pirimicarb concentration is raised from 0 to 50 mg/L (Fig. 4a).The pirimicarb binding to polymer caused the turn-on fluorescent intensity that is generally related to the transition from π to π* of the C=C bonds and benzene rings or from n to π* of C=O bonds or others, respectively.The fluorescence intensity was increased via www.nature.com/scientificreports/π-π interaction, and exhibits stable solid-state photo-luminescence (PL) to monitor the pirimicarb loading and releasing through the Förster resonance energy transfer with the pirimicarb like tramadol and doxorubicin [54][55][56] .The absorbance experiment showed that an interaction exists between pirimicarb and MIP polymer that turns on the absorbance behavior of MIP dependent on pirimicarb 57 .The fluorescence intensity of these samples at 395 nm are represented in Fig. 4c.The linear model of these data was calculated by Excel software (Eq.( 11)).As can be seen, the R 2 and sensitivity are equal to 0.9518 and 15.808 a.u.L/mg, respectively.The linear range can be measured by making a plot of analyte concentration versus fluorescence, and seeing at what concentration the data deviate from a straight line that is tangent to the low end of the concentration range.Therefore, the fluorescence intensity shows an excellent linear relationship with the pirimicarb concentration in 1 to 50 mg/L.
The results of measurement of real samples (Fig. 4b) are reported in Table 3.Also, the relative standard deviations (RSD) of this measurement were calculated.As can be seen, the RSD of all data is less than 2, so real samples can be measured with acceptable accuracy by this chemosensor.The LOD and LOQ of this sensor were calculated by the calibration line at low concentrations.The results show that the LOD and LOQ are equal to 1.79 mg/L and 5.43 mg/L of pirimicarb concentration, respectively.

Selectivity and binding sites
To examine the selectivity of NIP-CuFe 2 O 4 /GQDs and MIP-CuFe 2 O 4 /GQDs nanocomposites, competitive adsorption of pirimicarb was performed separately with each of Diazinon, Deltamethrin, and Chlorpyrifos pesticides, and results are reported in Table S1 (see SI).According to the results, the NIP nanocomposite has a virtually amount of distribution coefficient (K d ) for all intended pesticides.In contrast, the distribution coefficient of MIP nanocomposite to the adsorption of pirimicarb is considerably higher than that of the others.This reveals the high capability of selective pirimicarb adsorption by MIP nanocomposite, and a higher selectivity coefficient (K) for the MIP compared to the NIP nanocomposite.Besides, the relative selectivity coefficients of pirimicarb with each of the three other pesticides are roughly similar.In addition, this study was done for real samples and results of them are shown in Table S2.The relative selectivity coefficients of pirimicarb in real samples are similar to the results of laboratory samples.The Scatchard plot was also used to study the binding sites of MIP nanocomposite and compare it to NIP nanocomposite, and the results are shown in Fig. 5 and Table 4.As can be observed, the q max coefficient for MIP nanocomposite is higher than that of NIP nanocomposite, and (11) y = 15.808x+ 273.14  www.nature.com/scientificreports/ the imprinting factor value attained 2.48, which confirms the specific binding properties of pirimicarb pesticide.In Table 5, the comparison of the results of this study with other studies is reported.As can be seen, the chemosensor used in this study can compete with other similar chemical sensors.In addition, it has a much simpler structure than them.The detection limit value of the sensors reported in other studies is higher than the sensor developed in this study, but their working and linear ranges are also different.However, the maximum residue level of pirimicarb pesticide is equal to 3 ppm, which in this study; the developed sensor has the ability to detect much lower and higher amounts than.In addition, imprinting factor of MIP-CuFe 2 O 4 /GQDs is also suitable and comparable with other reported polymers.

Reusability
For investigating the reusability of MIP-CuFe 2 O 4 /GQDs, in the first step, the best solution was selected (13 different pH of Britton-Robinson buffer).The results have been shown in Fig. S5-a (see SI).As can be seen, the best solution is the Britton-Robinson buffer at pH = 1.Then, up to 15 cycles of pirimicarb adsorption by the MIP-CuFe 2 O 4 /GQDs were carried out, and after each step, the polymer was washed by Britton-Robinson buffer (pH = 1).After each cycle, the percentage of adsorption was calculated, and the results are shown in Fig. S5-b (see SI).After 15 cycles, the adsorption percentage of MIP-CuFe 2 O 4 /GQDs for pirimicarb is higher than 96%.Therefore, MIP-CuFe 2 O 4 /GQDs have good reusability for pirimicarb pesticide as an optical chemosensor.

Conclusion
In this study, the MIP-CuFe 2 O 4 /GQDs were synthesized and characterized by FT-IR, XRD, TEM, and SEM techniques.The results of this analysis showed the MIP-CuFe 2 O 4 /GQDs was correctly synthesized.Then, the ability of MIP-CuFe 2 O 4 /GQDs as an optical chemosensor was investigated.In the first step, the optimal condition for the adsorption of pirimicarb by this polymer was studied.The results showed the best conditions for maximum adsorption are pH = 7.54, initial pirimicarb concentration = 10.17 mg/L, polymer dosage = 840 mg/L, and contact time = 6.15 min.And in these conditions, the adsorption percentage is 99.92.The sensitivity, LOD, and LOQ for this optical sensor were calculated based on their fluorescence intensity, and are equal to 15.808 a.u.L/mg, 1.79 mg/L, and 5.43 mg/L, respectively.The selectivity and reusability of this sensor were studied, and the results demonstrated that the imprinting factor for MIP-CuFe 2 O 4 /GQDs is equal to 2.48.Reusability results of MIP-CuFe 2 O 4 /GQDs as a sensor shows that their adsorption percentage after 15 cycles is higher than 96%.Generally, according to the results of this study, the MIP-CuFe 2 O 4 /GQDs is the best chemosensor material for the detection of pirimicarb in an aqueous solution.
9) % of Desorption = amount of pirimicarb desorbed amount of pirimicarb adsorbed × 100 and MIP-CuFe 2 O 4 /GQDs nanocomposite, the distance between the crystal plates was determined by using the Bragg law, and the size of the crystals was calculated using the Debye-Scherrer relationship.According to the results, the crystal sizes for them are 26.5 and 13.2 nm, respectively, and the distance between the crystal plates is equal to 0.253 nm.The result of BET analysis is shown in Fig. 1c.As can be seen, the surface area of CuFe 2 O 4 NPs in BET isotherm is equal to 83.64 m 2 /g.The surface-to-volume ratio of CuFe 2 O 4 NPs is very large and more suitable compared with other nanoparticles such as iron nanoparticles.The BET surface area of Fe 2 O 3 NPs 44 and CuFe 2 O 4 NPs 45 approximately were equal to 30 m 2 /g and 63 m 2 /g, respectively.Vibrating sample magnetometer (VSM) analysis was used to investigate and determine the magnetic properties of CuFe 2 O 4 nanoparticles and MIP-CuFe 2 O 4 /GQDs nanocomposite (Fig. 1d).The results showed that the CuFe 2 O 4 nanoparticles and MIP-CuFe 2 O 4 /GQDs nanocomposite have hysteresis.The magnetic saturation values of CuFe 2 O 4 nanoparticles and MIP-CuFe 2 O 4 /GQDs nanocomposite are 33.20 emu g −1 and 8.31 emu g −1 , respectively.As can be seen, the magnetic saturation value of MIP-CuFe 2 O 4 /GQDs is very small compared to that of CuFe 2 O 4 nanoparticles, which indicates the coverage surface of the CuFe 2 O 4 with polymer, completely.Figure 2a-c show SEM images of CuFe 2 O 4 NPs, GQD, and MIP-CuFe 2 O 4 /GQDs, respectively.The SEM image of graphene quantum dot shows its fully layered structure.CuFe 2 O 4 nanoparticles have a spherical shape with an average diameter of about 25.97 nm (Fig. 2d).In the case of MIP-CuFe 2 O 4 /GQDs, the imprinted layer surrounded the CuFe 2 O 4 NPs, and the average particle size increased to about 255.87 nm (Fig. 2e).The TEM image for CuFe 2 O 4 NPs and MIP-CuFe 2 O 4 /GQDs nanocomposite are shown in Fig. 3.As can be seen in Fig. 3a, the CuFe 2 O 4 NPs are clear and Fig. 3b shows CuFe 2 O 4 NPs with a polymer layer on their surface.The SEM results are validated by the results of TEM.

Figure 4 .
Figure 4.The results of fluorescence intensity with pirimicarb concentration: (a) all of the curves for experimental samples, (b) all of the curves for real samples, (c) calibration curve.

Table 1 .
Independent variables and levels of RSM.

Table 3 .
Results of measurement of real samples and RSD.

Table 4 .
Results of Scatchard analysis.

Table 5 .
Comparison of the results with other studies.